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ABSTRACT 

Most type-Ic core-collapse supernovae (CCSNe) produce and neutron 
stars (NSs) or black holes (BHs). The dipole radiation of nascent NSs has usually 
been neglected in explaining supernovae (SNe) with peak absolute magnitude 
Mpeak in any band are > —19.5 mag, while the ^®Ni can be neglected in fitting 
most type-Ic superluminous supernovae (SLSNe Ic) whose Mpeak in any band 
are < —21 mag, since the luminosity from a magnetar (highly magnetized NS) 
can outshine that from a moderate amount of ^®Ni. For luminous SNe Ic with 
—21 < Mpeak ^ —19.5 mag, however, both contributions from ^®Ni and NSs 
cannot be neglected without serious modeling, since they are not SLSNe and the 
^®Ni mass could be up to ~ 0.5Mq. In this paper we propose a unihed model that 
contain contributions from both ^®Ni and a nascent NS. We select three luminous 
SNe Ic-BL, SN 2010ay, SN 2006nx, and SN 14475, and show that, if these SNe are 
powered by ®®Ni, the ratio of Mni to Mgj are unrealistic. Alternatively, we invoke 
the magnetar model and the hybrid (^®Ni -|- NS) model and find that they can fit 
the observations, indicating that our models are valid and necessary for luminous 
SNe Ic. Owing to the lack of late-time photometric data, we cannot break the 
parameter degeneracy and thus distinguish among the model parameters, but 
we can expect that future multi-epoch observations of luminous SNe can provide 
stringent constraints on ^®Ni yields and the parameters of putative magnetars. 
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Subject headings: stars: magnetars, - supernovae: general, - supernovae: individ¬ 
ual (SN 2010ay, SN 2006nx, SN 14475) 

1. Introduction 

It has long been believed that most aged massive (zero-age main-sequence mass Mzams ^ 
S.OMq) stars terminate their lives as core-collapse supernovae (CCSNe) (Woosley et ah 2002; 
Janka et ah 2007), which are classified into type IIP, IIL, Iln, Ilb, Ib and Ic according to their 
spectra and light curves (Filippenko 1997), leaving neutron stars (NSs) or black holes (BHs) 
at the center and producing a moderate amount of ^®Ni which is generally regarded as the 
dominant power source for most supernovae (SNe) (Colgate &: McKee 1969; Colgate et ah 
1980; Arnett 1982) Of all these subclasses, SNe Ic have attracted more and more atten¬ 
tions since a great number of these events have been discovered and confirmed in recent 
years that some SNe Ic with broad absorption line features (“broad-lined” or “BL”) have 
an accompanying gamma-ray burst (GRB) or X-ray hash (XRF) (Woosley &: Bloom 2006; 
Hjorth & Bloom 2012). Main models of central engines of GRBs associated with SNe are the 
“collapsar” model (Woosley 1993; MacFadyen & Woosley 1999) involving “BH -|- disk” sys¬ 
tems and the magnetar (highly magnetized neutron star) model (Usov 1992; Metzger et ah 
2007; Bucciantini et ah 2008; Metzger et ah 2011) proposing that the explosions may leave 
benind fast-rotating magnetars, being also regarded as the origin of the shallow decays and 
plateaus as well as rebrightenings in the multi-band afterglows of some GRBs (Dai & Lu 
1998a,b; Zhang & Meszaros 2001; Dai 2004; Dai &: Liu 2012). 

On the other hand, in the last two decades, optical—NIR and radio observations have 
revealed that not every SN Ic-BL associates with a GRB or an XRF (Soderberg et ah 2005; 
Drout et ah 2011). Many SNe Ic-BL, no matter whether they are associated a(n) GRB/XRF 
or not, have very high kinetic energy ^ 1-0 x 10®^ erg and have therefore been called 
“Hypernovae” (Iwamoto et al. 1998). Supposing that the optical—NIR emission is powered 
by radioactive ®®Ni decay, these SNe Ic need ~ 0.1 — O.5M0 of ®®Ni. All GRBs associated 
SNe Ic and most GRB-less SNe Ic are not very luminous, with peak absolute magnitude 
Mpeak > -19.5. 


^Since the observations for SN 1987A (Lundqvist et al. 2001) and SN 1998bw (Sollerman et al. 2002) have 
already revealed that the radioactive elements other than ®®Ni (e.g., ®^Ni, ^"^Ti and ^^Na, etc) constitute 
only a minor fraction of the radioactive masses and contribute dominant fluxes at very late times (>600 days 
and 1,200—1,400 day after explosion for SN 1987A and SN 1998bw, respectively), we only consider the 
contribution of ®®Ni, which is the most abundant nucleus resulting from explosive silicon burning in shock- 
heated silicon shells and the dominant energy source at early-time (e.g., < 500 days) of a SN. 
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Thanks to the unprecedented boom of targeted- and untargeted- sky survey programs, 
many superluminous SNe (SLSNe) whose Mpeak in any band are < —21 mag (Gal-Yam 
2012) have been found in the past decade, most of which cannot be explained by the 
widely adopted ®®Ni decay model (Quimby et ah 2011; Inserra et ah 2013; Nicholl et ah 
2013; McCrum et al. 2014; Nicholl et ah 2014), motivating researchers to consider alter¬ 
native energy-reservoir models. Currently, main models explaining the SLSNe are SN ejecta 
- circumstellar medium (GSM) interaction model (Chevalier 1982; Chevalier & Fransson 
1994; Chevalier & Irwin 2011; Ginzburg &: Balberg 2012) that has been employed to ex¬ 
plain many Type Iln luminous SNe (Chugai 1994; Zhang et al. 2012; Ofek et al. 2013) and 
superluminous SNe Iln (Smith & McCray 2007; Moriya et al. 2013) as well as some Type 
Ic SLSNe (Chatzopoulos et al. 2013; Nicholl et al. 2014), and the magnetar-powered SLSNe 
model (Kasen & Bildsten 2010; Woosley 2010) that has been used to explain many Type 
Ic SLSNe (Inserra et al. 2013; Nicholl et ah 2013; Howell et al. 2013; McCrum et al. 2014; 
Vreeswijk et al. 2014; Nicholl et al. 2014; Wang et al. 2015). In some cases, the magnetar- 
powered SLSNe model with the assumption of full energy trapping fails to £t the late-time 
light curves (Inserra et ah 2013; Nicholl et al. 2014). In Wang et al. (2015), we generalized 
the magnetar-powered SLSNe model by introducing the hard emission leakage and solved 
this problem, highlighting the importance of the leakage effect in this model. 

Thus, the power sources of most normal SNe and SLSNe have been attributed to ®®Ni- 
decay and the magnetar spin-down or ejecta-CSM interaction, respectively. Within this 
picture, the dipole radiation of nascent magnetars has usually been neglected in explaining 
the SNe with Mpeak ^ —19.5 (but see Maeda et al. 2007), while the ®®Ni decay energy 
have generally been ignored in fitting all Type Ic SLSNe (Mpeak ^ “21) because the peak- 
luminosity due to the energy injection from a newly born fast-rotating (initial period ~ 
millisecond) magnetar can outshine that from a moderate amount of ^®Ni (see Inserra et al. 
2013, for further analysis). 

Besides normal SNe and SLSNe mentioned above, however, there are some SNe with 
—21 < Mpeak ^ —19.5 mag, constituting a class of “gap-hller” events that bridge normal 
SNe and SLSNe, have been found by many telescopes. Most of these “gap-hller” SNe are 
explained by the ®®Ni decay model for SNe Ic as well as “Super-Chandrasekhar-Mass” SNe 
la, and the ejecta-CSM interaction model for SNe Iln. 

Among these “gap-hller” SNe with —21 < Mpeak —19.5 mag, luminous Type Ic 
SNe and their energy-reservoir mechanisms have not attracted enough attentions yet, their 
high peak-luminosities are simply attributed to the ®®Ni cascade decay without any detailed 
modeling. It has recently been demonstrated that the ®®Ni-decay model cannot be arbi¬ 
trarily used to explain all SNe Ic especially those very luminous ones (Inserra et al. 2013; 
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Nicholl et al. 2013) since the production of ®®Ni is ineffective in CCSNe and more ^®Ni mass 
needs more mass of the ejecta, which could theoretically result in broad light curves that 
usually conflict with observations. While the CSM-interaction model can predict a wide 
range of peak luminosities and therefore explain the normal, luminous, and superluminous 
SNe Iln as a whole category, ®®Ni-decay model cannot explain luminous SNe Ic. 

Here, we propose that for some luminous SNe Ic, there are other energy reservoirs 
that play a signihcant role in producing light curves. Since all hydrogen envelopes of pro¬ 
genitors of SNe Ic have been stripped and the spectra are lack of narrow and/or emission 
lines indicative of the interactions we can exclude the ionized hydrogen re-combination 
(Falk & Arnett 1977; Klein & Chevalier 1978; Popov 1993; Kasen & Woosley 2009) and do 
not consider the ejecta-CSM interaction process. Therefore, like the cases of SLSNe Ic, the 
nascent magnetar embedded in the center of the explosion is the most promising candidate 
for the additional power source (Ostriker & Gunn 1971). On the other hand, the production 
of 5®Ni is inevitable and cannot be neglected in modeling luminous SNe Ic because their 
peak luminosities are considerably lower than that of SLSNe and the contributions from 
a moderate amount of ^®Ni can be comparable with the contributions given by any other 
energy reservoirs. Furthermore, in principle, the contributions from ®®Ni and NSs cannot be 
directly neglected for CCSNe with a wide range of peak luminosities. Thus, it is necessary 
to consider a unihed model containing contributions from both ^®Ni and NSs. 

In this paper we construct the unihed model which contain the contributions from both 
the ®®Ni cascade decay energy and the NS rotational energy, and apply it to explain the 
light curves of some luminous SNe Ic-BL, i.e., SN 2010ay, SN 2006nx, and SN 14475. The 
remainder of this paper is organized as follows. In Section 2, we give a unihed semi-analytical 
model for Type Ic SNe. Based on this model, we ht the light curves of SN 2010ay, SN 2006nx, 
and SN 14475 in Section 3. Finally, some discussions and conclusions are presented in Section 
4. 


2. The unified semi-analytical model for SNe 

In this section, we construct a unihed semi-analytical model combining energy from a 
spin-down magnetar and an amount of ^®Ni to describe the SN luminosity evolution. In this 
unihed model, ®®Ni release high energy photons via the cascade decay chain ^®Ni—>^®Co^^®Fe, 
while the magnetar can inject its rotational energy into the SN ejecta as heat energy 


^the unique exception is SN 2010mb which cannot be explained by ^®Ni-powered model and has a signature 
of interactions (Ben-Ami et al. 2014) 
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(Ostriker & Gunn 1971). 

Based on Arnett (1982), taking into account the 7 -ray and X-ray leakage (e.g., Clocchiatti & Wheeler 
1997; Valenti et al. 2008; Chatzopoulos et ah 2009, 2012; Drout et ah 2013; Wang et al. 

2015), the luminosity is given by 

9 7 , ZflQtA pt f t '2 , 2Rot' 

L(t) = _e '"’’m J (x_ 

Tm ^ Jo 

X f ^ P{t')dt' erg s“\ (1) 

\VTra TmJ 


where Ro is the initial radius of the progenitor, which is very small compared to the radius 
of the ejecta. We take the limit Rq 0, then the above equation can be largely simplihed. 
With Equations (18), (19) and (22) of Arnett (1982), the effective light-curve timescale Tm 
can be written as ^ 


7 2KMej \ 
\^) ’ 


where k, is the optical opacity to optical photons (i.e., the Thomson electron scattering 
opacity), Mej and v are the mass and expansion speed of the ejecta, respectively, and c is 
the speed of light. /9 ~ 13.8 is a constant that accounts for the density distribution of the 
ejecta. P{t) is the power function. Here, v is the scale velocity (ugc) in Arnett (1982) and 
approximates to the photospheric expansion velocity Uph. Hereafter, we let v ~ Uph. 

The factors and (1 — in Equation (1) represent the 7 -ray leakage and 

trapping rate, respectively. T^(t) = At~‘^ is the optical depth to 7 -rays (Chatzopoulos et al. 
2009, 2012). If the SN ejecta has a uniform density distribution (Mej = ( 4 / 3 ) 7 rpi?^, 

Ek = (3/10)MejU^), A depends on (the opacity to 7 -rays), Mej and v as 


A 


SKjMgj 


4.75 X 10 


13 




0.1 cm^ g ^ 


X 


Me 


M 


o 


10 ® cm s“^ 


-2 


( 3 ) 


^The coefficient 2 in the equation has been adopted as ^ in some other papers (Chatzopoulos et al. 
2009, 2012, 2013; Wang et al. 2015), the latter originated from a typo in Equation (54) of Arnett (1982) that 
written 5/3 as 3/5, see also the footnote 1 in Wheeler et al. (2014). To get the same light curves, k or Mej 
(v) should be multiplied by 5/3 (3/5), or these three parameters simultaneously adjusted so that Tm can be 
invariant for an individual event. 
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2.1. The ^®Ni-decay energy model 

In the ®®Ni-decay energy model, the input power is 

g ^/tCo _ g ^/"T'Ji 

P{t) = Pm{t) = + ecoMNi— - - -erg s"\ (4) 

1 - Tm/Tco 

where eNi = 3.9 x 10^° erg s“^ g“^ is the energy generation rate per unit mass due to 
®®Ni decay (^®Ni^®®Co) (Cappellaro et ah 1997; Sutherland & Wheeler 1984), is the 
initial mass of ®®Ni, tni = 8.8 days is the e-folding time of the ^®Ni decay, eco = 6.8 x 10® 
erg s“^ g“^ is the energy generation rate due to ^®Co decay (^®Co —>■ ^®Fe) (Maeda et al. 
2003), Tco = 111.3 days is the e-folding time of the ^®Co decay. 


2.2. The magnetar spin-down energy model 


In the magnetar model, supposing that all the spin-down energy released can be con¬ 
verted into the heat energy of SN ejecta and the angle between the dipole magnetic helds 
and the magnetar’s spin axis is 45°, the power coming from the magnetar dipole radiation 
is (Ostriker & Gunn 1971) 

P{t) = PNs{t) = A -^ erg s-\ (5) 

Tns (1 V^ns) 

tns = 6/Nsc^/7?^i?Ns^o = 1-3 (5/10^'^ G)“^(Po/10 ms)^ yr is the spin-down timescale of the 
magnetar (Kasen & Bildsten 2010). is fhe rotational energy of the magnetar. 


PnS — (l/2)/NsflNS 

52 ^NS 


~ 2 X 10 


1.4 Mq \1 ms 


Pn 


-2 


Rns 

10 km 


erg. 


( 6 ) 


where Mns, Pns and Pq are the mass, radius and initial rotational period of the magnetar, 
respectively, while Ins = (2/5)MnsPns i® li^e moment of inertia of the magnetar whose 
canonical value is 10^^ g cnP (Woosley 2010). 


2.3. The hybrid (®®Ni -|- magnetar) model 

If an amount of ^®Ni is synthesized and a fast-rotating magnetar is left after a SN 
explosion, both contributions from these two power sources must be taken into account, this 
is the hybrid model. In this model, the luminosity of a SN is 


Lum{t) — Lj<ii(t) + LNs(i), 


(7) 







7 


where Lm{t) and -hNs(^) are luminosities supplied by ^®Ni and the magnetar, respectively. 

In this paper we propose the “unihed” model that contains the above three models, i.e., 
the ^®Ni model, the magnetar model, and the hybrid model. If the contribution from the 
magnetar or ®®Ni can be neglected, this model would be simplihed to the magnetar model 
or the ^®Ni-decay model, respectively. The parameter set for the unihed model is (Mej, v, 
Mni, B, Pq, When M^i = 0, the parameter set is (Mej, v, B, Pq, corresponding 
to the magnetar model; when the neutron star is non-rotational, the parameter set is (Mej, 
V, Mni, k^), corresponding to the ®®Ni-decay model; when the contributions from ®®Ni and 
the magnetar are both important, the model can be termed as the hybrid model. It should 
be noted that if the photospheric velocity Uph is not measured, then v ~ Uph is also a free 
parameter. 


3. Fits to the light curves of SN 2010ay, SN 2006nx, and SN 14475 

To verify this unihed model, we select three luminous SNe Ic-BL, SN 2010ay, SN 2006nx, 
and SN 14475. SN 2010ay was discovered by the Catalina Real-time Transient Survey (CRTS; 
Drake et ah 2009) with V-Band and R-Band peak absolute magnitudes. My ~ —19.4 mag 
and Mfi ^ —20.2 ± 0.2 mag, respectively. We use the R—band light curve as a proxy for the 
bolometric light curve of SN 2010ay. Although the former is not a completely reliable proxy 
for the latter in some cases (e.g.. Walker et ah 2014), Wheeler et ah (2014) pointed out that, 
for SN 1993J (Type Ilb), SN 1998bw (Ic-BL) and SN 2002ap (Ic-BL), the discrepancies 
between the R—band light curves and the quasi—bolometric light curves are rather small 
and can be neglected. 

SN 2006nx and SN 14475 were discovered by the SDSS-II (Taddia et ah 2015) with peak 
absolute bolometric magnitude ~ —20.36 mag and ~ —20.1 mag, respectively. Taddia et ah 
(2015) have constructed their bolometric light curves, so hereafter we adopt their data. 

The peak luminosities of SN 2010ay, SN 2006nx, and SN 14475 are ~ 4.2, 4.84, and 3.8 
times that of SN 1998bw which is a well-studied SN Ic-BL with peak absolute magnitude 
~ —18.65 and the ^®Ni mass ~ 0.43lg;o5 Mq (Mazzali et al. 2006), respectively. Hence, 
SN 2010ay, SN 2006nx, and SN 14475 are luminous SNe/HNe and brighter than all GRB- 
associated SNe/HNe conhrmed to date. 

The so-called Arnett law (Arnett 1979, 1982) reports that the peak luminosity of 
a SN purely powered by ®®Ni decay is proportional to the instantaneous energy deposi- 


tion, therefore is proportional to the initial ^®Ni mass which is widely invoked to infer 
the ^®Ni yields of some SNe (Contardo et ah 2000; Strolger et ah 2002; Candia et ah 2003; 
Stritzinger & Leibnndgnt 2005; Yuan et al. 2010; Cano et al. 2014). Comparing to the type 
Ic supernova SN 1998bw, the ®®Ni yields of SN 2010ay, SN 2006nx, and SN 14475 are es¬ 
timated to be l.SllgJJ Mq, 2.08lg;24 1.63lo;Jg respectively, when the Arnett 

law is applied. The value Mni ~ 1.811q;2} Mq for SN 2010ay is considerably larger 

than the value 0.9^q} Mq derived in Sanders et al. (2012) while the above values of Mjsa, 
2.08^0^4 SN 2006nx and 1.63lo;}g SN 14475 are roughly consistent with the 

values Mni = I.SOIq;}^ Mq and Mni = 1.27lo;o9 SNe derived in Taddia et al. 

(2015). 

The inferred ejecta masses of SN 2010ay, SN 2006nx, and SN 14475 are Mej ~ 4.7 Mq 
(Sanders et al. 2012), 7.52l^;74 Mq (Taddia et al. 2015), and 2.90l4;4g (Taddia et al. 
2015), respectively. Hence, the ratios of Mni to Mgj are ~ 0.4, ~ 0.2, and ~ 0.2 — 0.4 for 
SN 2010ay, SN 2006nx, and SN 14475, respectively. These values far exceed the values of 
all GRB-associated SNe, which are typically ~ 0.05 — 0.1 (Sanders et al. 2012) and larger 
than the upper limit (~ 0.2) for CCSNe (Umeda & Nomoto 2008). The above analysis 
suggests that there should be other energy sources aiding these three SNe to the high peak- 
luminosities. As mentioned above, power from ionized element re-combination and CSM- 
interaction can be safely neglected, we focus on the contributions of ®®Ni and the magnetar 
possibly leaved behind the stellar explosion. 

In this section, we use equations given in Section 2 to reproduce the semi-analytical light 
curves and £t the observations of SN 2010ay, SN 2006nx, and SN 14475. We assume that 
the hducial value of the optical opacity k is 0.07 (e.g., Taddia et al. 2015). The photospheric 
expansion velocity Uph of SN 2010ay, SN 2006nx, and SN 14475 are ~ 1.92 x 10*^ km s“^ 
(0.064 c) (Sanders et al. 2012), Ri 1.54 x lO'^ km s“^ (0.051 c), and ~ 1.87 x 10'^ km s“^ 
(0.062 c) (Taddia et al. 2015), respectively. Thus v is no longer a free parameter in our 
htting. 

In many light-curve modeling, has been assumed to have a hducial value 0.027 cm^ g“^ 


"‘Strictly speaking, the peak luminosity of a SN depends sensitively upon the opacity, initial ®®Ni mass, 
ejecta mass and kinetic energy. Only if the varieties other than ^®Ni mass are same for SNe or their effects on 
changing the peak luminosities of SNe can be canceled each other (e.g., higher opacity and lower ejecta mass, 
etc), then the peak luminosities of SNe are proportional to the initial ®®Ni mass. These conditions could 
be satisfied in many SNe, so the Arnett law can be used to infer the ^®Ni mass of SNe by comparing them 
to well-studies SNe (e.g., SN 1987A, etc) which have precise measurements of ^®Ni masses and bolometric 
light curves. In Section 3.1, we will use the ^®Ni-decay model to reproduce the light curves of SN 2010ay, 
SN 2006nx and SN 14475, and examine the ®®Ni masses derived by the Arnett law. 
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for the ^®Ni-powered SNe Ic (e.g., Cappellaro et al. 1997; Mazzali et ah 2000; Maeda et al. 
2003) and > 0.01 cm^ for the magnetar-powered SNe (see Kotera et ah 2013, Fig. 8), 
respectively. It is worth emphasizing that in the magnetar model, can vary between 
~ 0.01 and 0.2 cm^ g“^ {E^ > 10® eV) while kx can vary between ~ 0.2 and 10^ cm^ g“^ 
(10^ eV < Ex < 10® eV) (see Kotera et al. 2013, Fig. 8). Therefore, is also a free 
parameter in our htting. 


3.1. The ®®Ni-decay energy model 

The parameters for the ®®Ni-powered model are listed in Table 1 and the light curves 
reproduced by these sets of parameters are shown in Fig. 1 (Models Al and A2 for SN 2010ay, 
B1 and B2 for SN 2006nx, Cl and C2 for SN 14475). 

Using the empirical relation between the peak absolute magnitude and the mass of 
®®Ni derived by Drout et al. (2011), Sanders et al. (2012) inferred that the ®®Ni mass of 
SN 2010ay is O.OIq} Mq. However, it can be seen from Fig. 1 and Table 1 that 1.0 Mq 
of ®®Ni is inadequate to power the peak luminosity of SN 2010ay solely. To power the peak 
luminosity, ~ 2.0 Mq of ®®Ni must be required. The value 2.0 Mq is consistent with the 
value 1.811 o;2J Mq derived by the Arnett law. 

It can be seen from Fig. 1 and Table 1 that 2.0 Mq and 1.3 Mq of ®®Ni are required to 
power the peak luminosities of SN 2006nx and SN 14475, respectively. The values 2.0 Mq and 
1.3 Mq are consistent with the value I.SOIq;}^ Mq and 1.27^Q;Qg Mq derived by Taddia et al. 
(2015) as well as the values 2.08lo;24 Mq and 1.63lo;}g derived by the Arnett law. 

In this paper, we adopt a hducial value of optical opacity 0.07 cm^ g“^. In many other 
papers, however, the hducial values of opacity have also been assumed to be 0.06 cm^ g“^ 
(e.g., Valenti et al. 2011; Lyman et al. 2014), 0.08 cm^ g“^ (c-g-, Arnett 1982; Mazzali et al. 
2000), 0.10 cm^ g“^ (e.g., Nugent et al. 2011; Inserra et al. 2013; Wheeler et al. 2014) and 
0.2 cm^ g“^ (e.g., Kasen & Bildsten 2010; Nicholl et al. 2014). It is necessary to point out 
that more ®®Ni for the ®®Ni model or smaller Pq of the magnetar for the magnetar model is 
required to account for the peak luminosity of SN 2010ay, SN 2006nx, and SN 14475 when k 
is larger than 0.07 cm^ g“^ since larger k results in larger diffusion time and lower the peak 
luminosities. On the other hand, when k is smaller than 0.07 cm^ g“^, less ®®Ni or larger 
Pq are required for the same peak luminosity. For example. If we assume that SN 2010ay is 
powered by ®®Ni and k = 0.06 cm^ g“^, then ~ 1.8 Mq of ®®Ni must be synthesized, which 
is smaller than 2.0 Mq when k = 0.07 cnP g“^. 

To maintain the same peak luminosities and shapes of light curves, kMqJv must be 
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invariant. Becanse the photospheric velocities of SN 2010ay, SN 2006nx, and SN 14475 are 
fixed, so ^Mej =constant is required. Hence, smaller k requires larger Mej, as illustrated in 
Fig. 2. 


3.2. The upper limits of mass 


A consensus has been reached that iron—core-collapse SNe can synthesize a large amount 
of ^®Ni. However, could these 3 SN explosions produce ~ 1.3 Mq or ~ 2.0 Mq of ^®Ni? To 
get reasonable upper limits for the ^®Ni masses synthesized in the explosions of SN 2010ay, 
SN 2006nx, and SN 14475, we perform semi-quantitative estimates. 


The radius of the ^®Ni layer can be given as 


i?Ni — 


Eexp + E] 


bin 


jvraT'^ 


1/3 


~ 4.0 X 10^ 


Eem T F/bin \ 


1/3 


1051 gj-g J ^5 X 10® K 


-4/3 


km 


( 8 ) 


where E^^p, 77bin are the explosion energy and the binding energy of the progenitor, respec¬ 
tively, T is the temperature of the ejecta, a = 7.56566 x 10^^ erg cm“^ is the radiation 
energy density constant. Since Mni = 4/3pNi7ri?Ni) ^Ni is therefore approximately propor¬ 
tional to Eexp + F'bin — F'exp ^ which iu turn is approximately proportional to the kinetic 
energy E^, see, e.g., the bottom panel of Fig. 7 of Hamuy (2003) and the left panel of Fig. 
23 of Pejcha & Thompson (2014). The masses of the progenitors also influence the yields of 
^®Ni. These facts can also be seen in Fig. 7 and Table 3 in Umeda & Nomoto (2008). The 
typical masses of the ^®Ni synthesized during the explosive silicon burning are ~ 0.1 — O.5M0 
for CCSNe. 


We compare SN 2010ay, SN 2006nx, and SN 14475 with the prototype hypernova 
SN 1998bw. The kinetic energy of SN 2010ay (f» 1.1 x 10^^ erg, Sanders et al. 2012) and 
SN 14475 (Ri 4.7112^41® x 10^^ erg, Taddia et al. 2015) are about a factor of 3-5 smaller 
than that of SN 1998bw (Ri 3 — 5 x 10^^ erg, Iwamoto et al. 1998; Nakamura et al. 2001; 
Lyman et al. 2014), so we can infer that their ®®Ni yields must be < 0.5 Mq and ~ 
0.1 — 0.2 Mq which are signihcantly smaller than 1 Mq The kinetic energy of SN 2006nx 


^Pejcha & Thompson (2014) demonstrated that in most cases the absolute value of the progenitor binding 
energy i/bin must be less than the neutrino-driven wind energy E^nd (i.e., i/exp) and the total energy 
Etoti= Eexp + Ehin) is approximately equal to the wind energy i/wind, see their Fig. 12. 

® Although the kinetic energy of SN 14475 has rather significant uncertainty, its upper limit is only ~ 1.7 






11 


2.1601q 50^ X 10®^ erg, Taddia et al. 2015) has rather significant uncertainty, but its upper 
limit (Ri 5.91 x 10^^ erg) is slightly larger than that of SN 1998bw (r^ 5 x 10®^ erg). Hence, 
the ^®Ni yield of SN 2006nx must be ~ 0.5 Mq and smaller than 1 Mq. 

We can constrain the ^®Ni masses of SN 2010ay, SN 2006nx, and SN 14475 using another 
method. Umeda & Nomoto (2008) calculated the mass of the synthesized ^®Ni in CCSNe 
with progenitor Zero-Age Main-Sequence masses Mzams < IOOMq and found that the typical 
®®Ni mass is at most 20% of the ejecta mass. According to our light-curve modeling, the 
ejecta mass of SN 2010ay, SN 2006nx, and SN 14475 to be 6.5 Mq, 3.6 Mq, and 2.1 Mq, 
respectively, see Table 1. Combining these two results, < 1.3, 0.7, and 0.4 Mq of ®®Ni can 
be synthesized in the explosions, respectively. ^ If the masses of ^®Ni are ~ 2.0, 2.0, and 
1.3 Mq, the values of are ~ 0.31, 0.56, and 0.62 (see Table 1) which (far) exceed 

the upper limit (~ 0.2) of the mass ratio given by Umeda & Nomoto (2008) for CCSNe. 

These results indicate that the ^®Ni-decay model cannot account for the light curves 
of SN 2010ay, SN 2006nx, and SN 14475 and these three SNe are probably solely or partly 
powered by other energy sources. 


3.3. The magnetar spin-down energy model 

Encountered by the above difficulty, alternative energy-reservoir models must be seri¬ 
ously considered. According to their inferred ratio of Mni to Mgj, ~ 0.2, Sanders et al. (2012) 
suggested that SN ejecta-circumstellar medium interaction could also contribute to the high 
peak luminosity of SN 2010ay as well, yet they did not perform further investigation for this 
possibility. Owing to the lack of narrow and/or intermediate-width emission lines indicative 
of interactions between the SN ejecta and hydrogen- and helium-deficient CSM, we argue 
that ejecta-CSM interactions could hardly provide a reasonable interpretation for the excess 
luminosities of SN 2010ay, SN 2006nx, and SN 14475. Hence we attribute the luminosity 
excess instead to energy injection by the nascent magnetars and employ the magnetar model 
to explain the data for these SNe Ic. 

The parameters for the magnetar model are listed in Table 1 and the light curves repro- 


X10®^ erg, significant smaller than that of that of SN 1998bw, so we can expect that its ^®Ni mass should 
be smaller that of SN 1998bw. 

^It should be emphasized that these values are very rough and give rather loose upper limits. Nevertheless, 
real values of the masses of ®®Ni are difficult to be larger than them and the range of 0.1 —0.5 Mq is reasonable 
for these SNe. 
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duced by the sets of parameters are shown in Fig. 1 (Models A3 and A4 for SN 2010ay, B3 
and B4 for SN 2006nx, C3 and C4 for SN 14475). It can be seen from Fig. 1 that the mag- 
netar model can well £t the data. These resnlts indicate that this model can be responsible 
for the light curves of SN 2010ay, SN 2006nx, and SN 14475, i.e., SN 2010ay, SN 2006nx, 
and SN 14475 are probably powered by nascent millisecond (or ~ 10 ms) magnetars. 


3.4. The hybrid (®®Ni + magnetar) model 

When the ejecta masses of CCSNe are not very large, e.g., < 2 Mq, the ^®Ni synthesized 
are generally <0.1 Mq, which contribute a minor fraction of the luminosity of luminous SNe 
and can be neglected in the £t. Inserra et al. (2013) demonstrated that the magnetar models 
without ^®Ni and with ~ 0.1 Mq of ®®Ni are similar and therefore neglected the contribution 
from ^®Ni for their six SLSNe Ic. 

However, SN 2010ay, SN 2006nx, and SN 14475 are not SLSNe and their ejected ®®Ni 
might have a mass of ~ 0.1 — 0.5Mq, so the contribution from ^®Ni must be taken into 
account in our modeling. 

Since the exact values of Mni of these three SNe are unknown when the contribution 
from magnetars is considered, we plot light curves corresponding to a variety of masses of 
^®Ni. The typical values of adopted here are 0.1, 0.2, 0.5 and 1.0 Mq ®. The parameters 
for the hybrid (^®Ni + magnetar) model are also listed in Table 1 and a family of light curves 
reproduced by these sets of parameters are shown in Fig. 1 (Models A5—A8 for SN 2010ay, 
B5-B8 for SN 2006nx, C5-C8 for SN 14475). 

The ®®Ni mass with < 0.5 Mq, the energy released by the rapidly spinning magnetar is 
still the dominant energy source in powering the optical light curves of SN 2010ay, SN 2006nx, 
and SN 14475. When the ®®Ni mass is ~ 1 Mq, the energy released by ^®Ni is comparable 
with or exceed the energy released by the magnetar. If the ®®Ni masses are < 1 Mq and 
the magnetar parameters are adjusted, the light curves reproduced by the hybrid model are 
all in good agreement with observations, indicating that < 1 Mq of ^®Ni is admitted to be 
synthesized in their explosions. 

On the other hand, there are rather signihcant discrepancies in the late-time (> 50 days) 
light curves reproduced by different masses of ®®Ni. These discrepancies cannot be eliminated 
by adjusting the magnetar parameters since the decline rate of late-time light curves is 


^Although we have demonstrated in Section 3.2 that the ®®Ni mass of these three SNe can hardly be larger 
than 0.5 Mq, considering the possible uncertainties, we also adopt the value of 1.0 Mq in our modeling. 
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determined mainly by decay rate and the value of 

We cannot determine the precise masses of ^®Ni synthesized in the explosions of these 
SNe, because there is no precise late-time observation. Provided that there are some late¬ 
time data observed, it would be powerful enough to constrain the precise value of their ^®Ni 
mass. 

To compare these models, we calculate the values of y^/d.o.f for all these theoretical 
light curves see Table 1. Since the ®®Ni-powered models are disfavored in explaining 
these three SNe, as discussed above, we do not discuss them here even if they give smaller 
X^/d.o.f.. Therefore, we only consider the magnetar-powered models and hybrid models. It 
can be seen from the last column in Table 1 that Models A4 (magnetar-powered model with 
hard emission leakage), B7 (magnetar and 0.5 Mq of ^®Ni), and C5 (magnetar and 0.1 Mq 
of ^®Ni) have the smallest x^/d-o.f. Nevertheless, due to the absence of late-time data, we 
still cannot conclude that these models are most favorable ones. 


3.5. An analysis for the origin of the kinetic energy 

It is necessary to take into account the kinetic energy coming from the PdV work. The 
tapped rotational energy of a nascent magnetar (T^ns — 2 x 10 ^^ (-Po/1 ms)”^ erg) must be 
split into radiation energy (i^rad) and kinetic energy (T'K.mag)- The former heat the ejecta, 
while the latter accelerate the ejecta. When Pq ~ 9 ms, Pns — 3.0 x 10^° erg. Even if 
all of the rotational energy of the putative magnetar is converted into kinetic energy of the 
ejecta pK,mag ~ -Ens — 3.0 x 10®*^ erg, it is still far less than the total kinetic energy 
(Pk ^ 5 X 10^^ erg) of SN 2010ay, SN 2006nx, as well SN 14475 and can therefore be 
completely neglected. 

We now turn our attention to the neutrino-driven mechanism (Bethe & Wilson 1985). 
Hydrodynamic supernova simulations suggest that the neutrino-driven mechanism can pro¬ 
duce SN energies ^ 2 x 10^^ erg (Ugliano et al. 2012), far less than the observed kinetic 
energy of Pk ^ 5 x 10 ^^ erg, indicating that there must be another mechanism accounting 
for the additional energy of > 3 x 10^^ erg. A popular scenario to explain energetic SNe 


®Since Taddia et al. (2015) did not provide the observational errors of the data of SN 2006nx and SN 14475, 
we adopt a hducial value for the error, ^ 1.3 x 10^^ erg s“^, which is the smallest error of the SN 2010ay 
data, the difference between the real values of errors and our adopted values should not change the y^/d.o.f 
significantly. 

^°The calculations performed by Woosley (2010) have shown that i?k,mag — 0.4 when Pq = 4.5 ms, 
B = lx 1014 G. 
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with kinetic energies > 10^^ erg is the “collapsar” model involving a “BH-disk” system, in 
which the “BH-disk” system generates bipolar jet/outflow while the “disk wind” ensures that 
the progenitor star explodes and synthesizes enough amount of ^®Ni. However, our fittings 
have demonstrated that the newly formed magnetar did not collapse to a black hole at least 
several days after the explosion. Hence, the huge kinetic energy might stem from some mag¬ 
netohydrodynamic processes (e.g., magnetic buoyancy, magnetic pressure, hoop stresses, etc) 
related to the proto-NS themselves (Wheeler et ah 2000) after which the magnetars inject 
their rotational energy to heat the ejecta and generate the light curves. 


4. Discussion and Conclusions 

Most type Ic CCSNe produce ^®Ni which release energy to the ejecta via cascade decay 
and central NSs which convert a portion of their rotational energy into heating energy of 
the ejecta via the magnetic dipole radiation. For almost all normal CCSNe, ~ 0.1—0. 7 Mo 
of ^®Ni are adequate to power the peak and post-peak decline of observational light curves 
and the contribution of NSs can be neglected In contrast, for SLSNe Ic which are lack of 
evidence of interaction and cannot be explained by ^®Ni solely, millisecond NSs (magnetars) 
can supply almost all input power and the contribution of ^®Ni can be neglected. 

For luminous SNe Ic with —19.5 > Mpeak ^ ~21 mag in any band, however, both 
contributions from ®®Ni and NSs cannot be neglected without detailed modeling since the 
^®Ni permitted by completely explosive burning are usually inadequate to power the high 
peak luminosities while the production of ~ 0.1— ~ 0.5Mq of ^®Ni is reasonable and may 
contribute a significant portion of the total luminosity. Therefore, these “gap-filler” events 
which bridge normal SNe and SLSNe must be explained by a unified model which contain 
contributions from ^®Ni and NSs. 

To illustrate the necessity of constructing the unified model, we select and apply this 
model to three luminous SNe Ic-BL, SN 2010ay, SN 2006nx, and SN 14475. We tested 
the possibility that the luminosity evolutions of SN 2010ay, SN 2006nx, and SN 14475 are 
solely governed by radioactive ®®Ni cascade decay using the semi-analytical method and 
demonstrated that too large amounts, roughly 2.0, 2.0, and 1.3 Mq of ^®Ni, are needed 
to power their peak luminosities, respectively. These results disfavor the ®®Ni model in 


the explosion leaves a BH or nothing (corresponding to the so-called “pair instability SNe” 
(Heger & Woosley 2002) rather than CCSNe), then no NS contributes luminosity to the SN and the SN 
is solely powered by ®®Ni; if the explosion produces a transitional NS lasting about several hours or days the 
contribution of the NS to the light curve of the SN can also be neglected. 
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explaining the light curves of SN 2010ay, SN 2006nx, and SN 14475, since the ratios of Mni 
to Mej are 0.31, 0.56, and 0.62, respectively, signihcantly exceeding the upper limit (~ 0.2) 
calculated for CCSNe (Umeda & Nomoto 2008). 

Then we alternatively invoke the magnetar model to reproduce the light curves of 
SN 2010ay, SN 2006nx, and SN 14475 and show that the light curves reproduced can be 
well consistent with the observations, suggesting that SN 2010ay, SN 2006nx, and SN 14475 
are probably powered by spin-down energy injection from newly born millisecond magnetars. 
It seems that the magnetar model can account for the light curves of these SNe and ^®Ni 
is not necessary to be introduced. However, a moderate amount of ^®Ni, ~ 0.1 — O.5M0, 
can be expected in completely explosive burning of silicon shell and could give rather energy 
to those not very luminous SNe. So ^®Ni cannot be neglected when the SNe are not very 
luminous. Thus we use the hybrid model to £t these three SNe. 

We set = 0.1, 0.2, 0.5 and 1.0 Mq and hnd that < 0.2 Mq of ®®Ni can not sig¬ 
nihcantly inhuence the theoretical light curves while > 0.5 Mq of ®®Ni can effectively alters 
the shapes of the light curves, leading us to adjust the magnetar parameters to fit the ob¬ 
servations. Therefore, all these models (^®Ni, magnetar, and hybrid), with different amount 
of ^®Ni (from 0 to 2 or 1.3 Mq), £t the light curves equally well and can be unihed into a 
unihed model. Among these models, the ®®Ni model needs huge, unrealistic amounts of ®®Ni, 
resulting in over-high ratios of M^j to Mgj. In contrast, the magnetar model needs no ^®Ni 
and conhict with the basic theory of nuclear synthesis. The hybrid model is reasonable in 
explaining these luminous SNe. 

These facts suggest that these three luminous SNe Ic are all powered by both ®®Ni and 
magnetars, indicating that the magnetars contribute a large fraction of but not all power in 
driving their light curves. So we can conclude that, when Mpeak (in any band) is < —21 mag, 
the contribution of ®®Ni can be neglected in light-curve modeling; when —21 < Mpeak ^ 
— 19.5 mag, the contribution from the magnetar cannot be omitted in modeling while the 
contribution from ^®Ni must be considered seriously; but when Mpeak ^ —19.5 mag, the 
contribution of the magnetar can generally be neglected (but see Maeda et ah 2007). 

It should be pointed out that the “hducial cutoffs” —19.5 mag and —21 mag that divide 
the normal SNe, luminous SNe and SLSNe are artihcially given and do not have unambiguous 
physical meanings. Although no consistent picture of these SNe Ic with rather different peak 
luminosities has emerged so far, it seems that they have similar physical nature and similar 
power sources, e.g., ^®Ni and newly born NSs (magnetars), while the former result in a relative 
narrow range of peak-luminosities owing to the existence of the upper limit of ^®Ni yields 
and thus the peak luminosities, the latter might produce a wide range of peak-luminosities 
since the period of newly born NSs can vary from several milliseconds to several seconds. 
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and magnetic field strength of them can be up to 10^^ G. These two power sources could 
result in a continuous sequence of peak-luminosities, covering normal, luminous, and very 
luminous SNe Ic and unifying most of them into a whole category in which normal SNe 
are powered by slow-spinning NSs and ^®Ni, while luminous SNe and SLSNe are mainly 
powered by fast-spinning (millisecond) magnetar and ^®Ni. The similarity between late-time 
spectra of SLSNe Ic and spectra of normal SNe Ic-BL (Pastorello et ah 2010; Gal-Yam 2012; 
Inserra et ah 2013) also supports this unihed picture. 

Furthermore, we suggest that many SNe Ic with Mpeak ^ —19.5 mag might be partly 
powered by a nascent magnetar, i.e., the very young SN remnants of the explosions may 
harbor magnetars. The discrepancy in brightness between luminous SNe Ic and SLSNe Ic 
stem mainly from the initial spin periods of the magnetars, the former have Pq ~ 7 — 15 ms, 
while the latter have Pq ~ 1 — 6 ms (Inserra et ah 2013; Nicholl et ah 2014). For luminous 
SNe Ic, power from magnetars with Pq ~ 10 ms can become comparable to or exceed that 
from ^®Ni, and vice versa; for SLSNe Ic, millisecond magnetars usually overwhelm ^®Ni in 
powering them. 

For the sake of completeness, both these two power sources must be taken into account in 
light-curve modeling. Nevertheless, we still emphasize that this unihed scenario is especially 
important for the gap-hller SNe Ic since they, as aforementioned, are hardly powered solely 
by ®®Ni but ^®Ni can contribute a non-negligible portion of the luminosity and must be 
considered in modeling. Modeling for these luminous SNe can also help us to understand 
how normal SNe Ic relate to SLSNe Ic. Besides, even if a SN or superluminous SN that can 
be solely explained by the magnetar model, Fe line in its spectrum, if observed, must be 
explained by introducing a moderate amount of ^®Ni, because the magnetar model cannot 
explain Fe line (Kasen & Bildsten 2010). This is another advantage of the unihed model. 
In the recent excellent review for SLSNe, Gal-Yam (2012) already noticed these “gap-hller” 
SNe and reckoned that they are likely to be intermediate events between radioactivity- 
powered SNe/SLSNe Ic and the SLSNe Ic powered by some other processes. Our analysis has 
conhrmed this conjecture and furthermore demonstrated that these “intermediate events” 
are likely to have the same power sources and the diherence between them and their lower- 
luminosity and higher-luminosity cousins are mainly due to the physical properties of the 
nascent neutron stars embedded in their debris. 

While this unihed energy-reservoir model seems reasonable for SN 20I0ay, SN 2006nx, 
and SN 14475, the precise value of ^®Ni yields cannot be determined, indicating that the 
traditional method that suppose all power given by ^®Ni cascade decay is not always reliable 
in inferring the ^®Ni yields of GGSNe, especially for those luminous and very luminous ones. 
To discriminate these models and roughly determine the ^®Ni mass as well as the parameters 
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of the putative NSs, late-time (> 100 — 200 d) data are required. Unfortunately, these 
three SNe are all lacking late-time observations, especially for SN 14475 whose luminosity 
measurements in the period between 25 and 35 days after the peak are rather uncertain. 

Spectral analysis would be of outstanding importance in determining the precise values 
of ®®Ni and Mej as well as the parameters of the NS parameters. In practice, detection and 
analysis of nebular emission from ^®Fe resulted from the ^®Ni cascade decay and other heavy 
elements are necessary for precise measurement of ^®Ni yields. The ejecta mass Mej can 
signihcantly influence the rise time, shape, and peak luminosity as well as the ratio of Mni 
to Mej of a SN, so it is also a very important quantity in modeling. Due to the uncertainty of 
optical opacity k, Mej cannot be precisely determined using solely the light-curve modeling 
(see Fig. 2). To get more precise value of the ejecta mass or at least a rigorous upper/lower 
limit, as performed by Gal-Yam et ah (2009) for SN 2007bi, nebular modeling are required. 

The above procedures need very precise observations and spectral analysis lasting several 
hundred days after the detection of the hrst light. Unfortunately, no late-time photometric 
data and spectra have been observed for these three SNe, especially for SN 2006nx and 
SN 14475. We cannot get enough information to determine the precise masses of ^®Ni and 
the real values of magnetar parameters. Hence, we cannot distinguish among these model 
parameters and therefore to have a more precise scenario. We can expect that future high 
cadence multi-epoch UV—optical—NIR observations of the luminous SNe can pose stringent 
constraints on the parameters of the nascent NSs and ^®Ni yields. 
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Table 1: Parameters for unified (^®Ni-decay, magnetar spin-down, and hybrid (^®Ni -|- mag- 
netar)) modeling 



Me, 

B 

P 

V 


Ky 

Mni 

MNi/Mej 

xVd-o.f 


(A4) 

o 

Q 

(ms) 

(c) 

(cm^ g“^) 

(cm^ g“^) 

(Me) 



SN 2010ay 

Model Al “ 

6.5 

0 

oo 

0.064 

0.07 

oo 

2.0 

0.308 

5.41 

Model A2 “ 

6.5 

0 

oo 

0.064 

0.07 

0.014 

2.0 

0.308 

1.05 

Model A3 ^ 

6.5 

4.5 

00 

bo 

0.064 

0.07 

oo 

0 

0 

2.36 

Model A4 ^ 

6.5 

4.5 

00 

bo 

0.064 

0.07 

0.021 

0 

0 

1.34 

Model A5 

6.5 

4.5 

9.0 

0.064 

0.07 

0.021 

0.1 

0.015 

2.00 

Model A6 

6.5 

4.5 

9.3 

0.064 

0.07 

0.021 

0.2 

0.031 

2.05 

Model A7 

6.5 

5.2 

10.3 

0.064 

0.07 

0.021 

0.5 

0.077 

2.11 

Model A8 

6.5 

7.2 

12.5 

0.064 

0.07 

0.021 

1.0 

0.154 

2.12 

SN 2006nx 

Model B1 “ 

3.6 

0 

oo 

0.051 

0.07 

oo 

2.0 

0.556 

1.24 

Model B2 “ 

3.6 

0 

oo 

0.051 

0.07 

0.037 

2.0 

0.556 

1.22 

Model B3 ^ 

3.6 

4.5 

00 

bo 

0.051 

0.07 

oo 

0 

0 

1.58 

Model B4 ^ 

3.6 

4.5 

00 

bo 

0.051 

0.07 

0.037 

0 

0 

1.56 

Model B5 

3.6 

4.5 

9.0 

0.051 

0.07 

oo 

0.1 

0.028 

1.60 

Model B6 

3.6 

4.5 

9.2 

0.051 

0.07 

oo 

0.2 

0.056 

1.44 

Model B7 

3.6 

4.5 

9.8 

0.051 

0.07 

oo 

0.5 

0.139 

1.11 

Model B8 

3.6 

4.8 

11.9 

0.051 

0.07 

oo 

1.0 

0.278 

1.60 

SN 14475 

Model Cl “ 

2.1 

0 

oo 

0.062 

0.07 

oo 

1.3 

0.619 

4.10 

Model C2 “ 

2.1 

0 

oo 

0.062 

0.07 

0.063 

1.3 

0.619 

2.78 

Model C3 

2.1 

4.8 

10.4 

0.062 

0.07 

oo 

0 

0 

3.03 

Model C4 

2.1 

4.8 

10.4 

0.062 

0.07 

0.133 

0 

0 

3.30 

Model C5 

2.1 

4.8 

10.7 

0.062 

0.07 

0.133 

0.1 

0.048 

2.97 

Model C6 " 

2.1 

4.8 

10.9 

0.062 

0.07 

0.133 

0.2 

0.095 

3.53 

Model C7 ^ 

2.1 

4.8 

12.6 

0.062 

0.07 

0.133 

0.5 

0.238 

3.54 

Model C8 " 

2.1 

4.8 

15.2 

0.062 

0.07 

0.063 

1.0 

0.476 

4.03 


a rjpg niodel. 

^ The magnetar model. 

The hybrid (®®Ni -|- magnetar) model. 
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Fig. 1.— Model fits using the ®®Ni-decay models, the magnetar models as well as the hybrid 
(56Ni_decay + magnetar) models for SN 2010ay {top), SN 2006nx (middle), and SN 14475 
(bottom). Left panels: -10—100 days after the explosion; Right panels: -10—400 days after 
the explosion. Data and upper limits of SN 2010ay are obtained from Sanders et ah (2012), 
data of SN 2006nx and SN 14475 are obtained from Taddia et ah (2015). The horizontal 
axis represents the time since the explosion in the rest frame. Parameters for Models A1-A8, 
B1-B8, and C1-C8 are shown in Table 1. k = 0.07 cm^ g“^ is adopted. In all panels, we add 
two light curves reproduced by the ^®Ni-decay models with 1 Mq and 0.7 Mq, without and 
with 7 -ray leakage. It can be seen from these light curves that 1 Mq or 0.7 Mq of ®®Ni is 
inadequate to power these SNe. 
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K(cmV^) 


Fig. 2.— The values of k vs. the values of Mej for SN 2010ay, SN 2006nx, and SN 14475. 
The values of k vary from 0.05 cm^ to 0.20 cm^ g“^, resulting in rather different values 
of Mej required. Big Circle, Star, and Square correspond to SN 2010ay, SN 2006nx, and 
SN 14475 with k = 0.07 cm^ g“^, respectively. 





























































